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1. Introduction 

Chiral Perturbation Theory (CHPT) [1] is the modern way of exploiting Chiral Sym- 
metry constraints on strongly interacting processes. The mesonic two-flavour sector 
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was treated in an extensive work to next-to-leading order [2] . The physical processes 
and amplitudes treated in that reference were F-^, M^, the scalar and vector form 
factors of the pion, the pion radiative decay (tt —>■ iwj) and tttt scattering. Due to 
recent and planned experimental improvements, more accurate theoretical analyses 
are needed. The next-to-next-to-leading order of F^, were calculated in [3, 4], 
the pion radiative decay in [5] and tttt scattering in [4, 6]. In addition, the pion 
polarizabilities and the production of pion pairs in gamma-gamma collisions are also 
known to this order. For the latter, the neutral process was studied in [7] and the 
charged case in [3]. In this paper we present the full 0{p^) calculation of the scalar 
and vector pion form factors, thereby completing the calculation to the next order 
of the processes considered in [2]. 

In addition to the above calculations there are those where the amplitude is 
calculated using dispersive methods. This method allows for the full calculation 
up to the subtraction constants. These have then to be determined by comparison 
with experiment. The disadvantage of this approach is that one cannot do a simple 
comparison to existing models of low-energy constants appearing in CHPT, since the 
"chiral logarithm" parts of the subtraction constants cannot be fully determined^. 
As a consequence, we can neither vary the quark mass to compare with e.g. lattice 
calculations. The dispersive calculation for the form factors was done in [8], (see 
[9] for the analogous nn calculation) and an analytical expression for the dispersive 
integrals in the vector form factor was found in [10]. The one-loop formula and the 
latter partial two-loop results for the vector form factor have been used in the inverse 
amplitude method to fit over a larger kinematic range [11]. A different resummation 
scheme using constraints from Vector Meson Dominance (VMD), 1/A^c and unitarity, 
works as well in extending the range of validity [12]. For a recent review of the form 
factors in general see Ref. [13]. 

We have performed the calculation in two different ways. We have used the 
master formula approach as described in [4] and we have also directly computed 
all the relevant Feynman diagrams. Several nontrivial checks on the calculation of 
both form factors were done. The first is the absence of nonlocal divergences, i.e. the 
dependence on quark masses and external kinematical variables of the divergent parts 
must be analytic. Another powerful check results from the value of the form factors 
at zero momentum transfer. For the vector form factor this must be one because 
of gauge invariance and, for the scalar form factor, it is related to the derivative of 
the pion mass w.r.t. the quark mass. The latter relation follows from the Feynman- 
Hellman theorem. 

For both form factors we have abundant experimental information to compare 
with. In the vector case this information consists of direct measurements of the form 

^Sometimes the requirement that the limit gives a finite result allows to determine 

these terms, see [8] for an example. 
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factor, both in the timchke and spacclikc region. In the scalar case, since there is 
no microscopic scalar probe available, it is impossible to directly measure the form 
factor. On the other hand, it can be shown that the experimental information on 
the scalar, I — tttt phase shifts, inserted in a dispersive representation for the 
form factor, allows one to completely reconstruct the energy dependence of the form 
factor, modulo a multiplicative overall factor. Therefore in both cases we can make a 
detailed comparison of the CHPT two-loop expressions to the available experimental 
information. The fact that the form factors have a simple kinematical structure, 
makes even the two-loop representation rather simple and easy to manipulate. For 
this reason the comparison to the experimental information is particularly instructive. 



The contributions at order can be split into two different pieces: a dispersive 
contribution and a polynomial part. The numerical contribution of the dispersive 
part has been already analyzed in Ref . [8] . Inside the polynomial part we have again 
two types of terms: chiral logarithms and new 0{p^) low energy constants (LEG). 
The splitting between these two types of terms is arbitrary and depends on the 
renormalization scale /i: on the other hand we have learned from the experience at 
order that the p mass scale is a sensible choice for understanding the physical origin 
of these two types of terms. As we will show in our analysis, the same choice seems 
to be still sensible at order p^, since we will be able to understand at least the order 
of magnitude and sign of the new LEC's with an estimate based on the resonance 
saturation hypothesis. We will also confirm at order p^ the Vector Meson Dominance 
hypothesis, showing the different importance of the resonance contribution in the 
scalar and in the vector channel. Concerning the vector form factor, we stress that 
the direct comparison of CHPT to the data allows a rehable, model-independent 
extraction of the value of the charge radius. 



The paper is organized as follows. In Sect. 2 we define the form factors and 
the notation used. In Sect. 3 we describe the calculation and the checks performed, 
and give the analytical results for the form-factors and the associated radii. Sect. 4 
describes various estimates of the 0{p^) constants appearing in the calculation. In 
Sect. 5 we describe the other numerical input used and fit accurately to the available 
data for both form factors. Here we also describe the variation of the fits with sev- 
eral different assumptions and provide numerical results for all quantities, including 
F^/F, M^/M^ and an improved estimate of the low energy hadronic vacuum polar- 
ization contribution to the muon anomalous magnetic moment. In this section we 
also discuss the Omnes representation. We finally recapitulate our main conclusions 
in Sect. 6. 
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2. Definitions and Notation 

2.1 Form factors 

The scalar and vector form factors are defined respectively by 

{7r\p2)\uu + dd\7i^{p,)) = S'^Fsis^) , 
{Ap2)\1 {ui^u - d-f^d) |7r^(pi)) = ie'^^{pi^ + P2i.)Fv{s^) , (2.1) 

where s^^ = (p2 ~ Pi)^- The scalar form factor is defined through the isosp in-zero 
scalar source. The isospin-one scalar form factor can be defined analogously but 
it only starts at 0{p^). The vector form factor is also an isovector, and what we 
calculate here is its Iz — component. Similar definitions exist for the other isospin 
components. In what follows we take the isospin limit i.e. = = m, so that 
the other isospin-one components of the vector form factors are identical to the one 
we analyze here. 

2.2 Chiral Perturbation Theory 

Quantum chromodynamics with two fiavours in the chiral limit has an SU{2)l x 
SU{2)ji = 0(4) symmetry which is spontaneously broken down to SU{2)v = 0(3). 
At low energies, the three resulting Goldstone Bosons are the only relevant degrees of 
freedom and their interactions are strongly constrained by the underlying symmetry. 
The three Goldstone Bosons can be identified with the pions and the way of extracting 
the consequences of the chiral symmetry is Chiral Perturbation Theory. We use the 
non-linear sigma model or O (4) parametrization with the external field formalism of 
Ref. [2]. 

To lowest order in the low energy expansion, O(p^), processes are described by 
the tree-level diagrams of the lagrangian 

A = ^ V^C/t V^C/ + F^x'U) , (2.2) 
with the covariant derivative defined by 

V^[/° = a^C/° + 4(a;)C/\ 

V^W = d^W + e'^'vl{x)U' - al{x)U\ (2.3) 
where U{x) is an 0(4) four-component vector 

Here v^{x) and a^{x) are the external vector and axial-vector sources respectively and 
X — 2S(s°,p*) contains the isospin-zero part of the scalar source and the isospin-one 
part of the pseudoscalar source. 
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To O(p^) the amplitudes are given by one-loop graphs with vertices from £2 
and tree-level graphs containing vertices from L2 and one vertex from the O(p^) 
Lagrangian given by 



(2.5) 



where the strength tensor F^,^ is defined by 



(2.6) 



and X = 2B(p^,s'). 

The C(p^) contributions contain pure two-loop diagrams with vertices from £2, 
one-loop diagrams with vertices from £2 and one vertex from £4, and tree level 
diagrams with vertices from £.2 and either two vertices from £4 or one vertex from 
the 0{p^) Lagrangian. The latter Lagrangian, Cq, can be found in [14]. 

3. The calculation 

The calculation has been made using two different methods: in one case, we have 
calculated directly all the relevant Feynman diagrams with full generality. As an 
alternative method we use the master equation approach [4], which corresponds to 
recognizing that a large part of the graphs comes together such that they are one- 
loop graphs with one of the vertices given either by the one-loop scalar or vector 
form factor or by the tttt scattering amplitude with the pion legs off-shell. Some 
two-loop diagrams with non-overlapping loops have then to be evaluated separately 
in order to obtain the correct normalization. The integrals have been evaluated using 
the methods of [15], see also [3, 4, 5]. The subtraction procedure we used is a version 
of the modified minimal subtraction (MS) as described in [4]. 

Both methods have been used independently and yielded the same result. They 
also satisfy the requirements of gauge invariance, i.e. -Fy(O) = 1, and of the Feynman- 
Hellman theorem, i.e. -Fs(O) = dM^/drh. Both constraints couple quite a few 
diagrams in each process providing thus a good check on the calculations. In addition, 
at two-loop order there are nonlocal divergences diagram by diagram, which also 
cancel in the sum. The "double chiral logs" also satisfy the constraints imposed by 
renormahzation, meaning that the terms of the type (log(M|//x^))^ can always be 
cast inside the ki quantities defined in Sect. 3.1- see [16] for further explanation and 
references. 
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3.1 and 

For the sake of completeness we write down here the pion decay constant and the pion 
mass. This also serves to introduce the notation we use for the other quantities of 
interest. The difference with the expressions used in [3, 4] is that we have subtracted 
the infinities using MS scheme and rewritten the O(p^) part in terms of the physical 
mass and decay constant of the pion (M^ and respectively). 



F 



1 + X2ill - L) + xl 
7 



N 



12 



13 1 
192 iV2 



+k2- 21111 + 2{llf - -h + 



+ 0{xl) , 



(3.1) 



and 



r2 
M2 



-| = l + x,(2ZS + -L) + 



1(/^ + 2/-^l) + 



— h - 2h - + 4/3^4 - Ih + + r'j^ 



163 1 
'96 

+ 0{xl 



(3.2) 



The constants Vp and r]^ denote the contributions from the 0{p^) lagrangian after 
modified minimal subtraction. 

We have defined the following quantities 



N = 


167^^ 




X2 = 


Ml 

TT 




L = 




Ml 


ki = 


(4^[ - 






2Bm, 





(3.3) 



being the lowest order pion mass and F the pion decay constant in the chiral limit. 
The are the finite part of the coupling constants in £4 after the MS subtraction, 
and their values depend on the renormalization scale /j, as iJ,^{dll /dn^) = —{'^i/2N). 
The 7i were calculated in [2] and are given by 



7l = g; 72 = g; 73 



"2' 



1 1 

2; 75 = -7; 76 = 77 

6 3 



0. 



(3.4) 



Later we will also follow common use, and discuss the 0{p^) LEC's in terms of their 
scale-invariant combinations, the Zj's, defined as 



/I 



2N 



ik + NL) . 



(3.5) 
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3.2 Scalar Form Factor 



Wc start with the expression of the scalar form factor to two loops evaluated with 
the methods described above. 



Fs{s) = Fs{0) [1 + X2 {^{2s - 1) J» + s{ll -L-j^: 
+0{xl) . 



X2 



(3.6) 



As has been already mentioned, the two-loop contribution has been split into two 
parts: the polynomial piece of the amplitude reads 



11 
12 

32 



^ArU728 9 ^ 9^ " 



36 



1817 1 , 
1296 + 



31, 17, , 1 
_k, + -k,-k,^-^ 864 



9 



9 



9 



20 1 
'8TiV2 



(3.7) 



where 

s = sjMl ■ 

while the dispersive piece can be cast in the following form 



Uf\s)^j{s) 



-II (lls^ - 40s + 44) + \ll (75^ - 20s + 28) 
3 3 



+5/^ + Z[s^ + -s-l]+ —L -43s^ + 53s - 



18 



1 /29 2 61 391 

H — — s s H 

N V12 9 36 

--3(^)(r-i)- 



119 

~Y 
43 
36' 



4 1 

3^+4 



(3.8) 



The integral functions J, Ki, K2, and K4 are defined in [4] and we reproduce 
them here for the sake of completeness 



fj\ / 

K2 
\KsJ 



z -4N\ ( 
0^0 
8 



h 

V-(27V2)-V 



(3.9) 



and 



K4 = — { -Ki + -K3 H J + 

s^ 1 2 ^ 3 ^ 12N^ 



(3.10) 
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where 



h{s) 



:ln 



1 - 



(3.11) 



v^+1 ' ' s 

The functions s~^J and s'^Ki are analytic in the complex s-plane (cut along the 
positive real axis for s > 4), and they vanish as |s| tends to infinity. Their real and 
imaginary parts are continuous at s = 4. 

We use the form factor at zero momentum transfer given by 



Fs{0)=2bI^1 + x, (4/^ + L + ^) 

- f-11/: - 21', - 31; + 11- —l] + — At 



+xl 



21 21 1 



+ 0{xl) , (3.12) 



as a check of the previous result, Eq. (3.6). One can see that this last expression 
can also be derived using Eq. (3.2) and the Feynman-Hellman theorem. It agrees 
exactly and leads to the relation between the 0{p^) constants r^^ = 3r]^. 

We can now expand the form factor for s -C 1 (stt -C 4M^), obtaining the 
expression 

Fs{s) = Fs{0) (l + lir^ss + c^s' + 0(.^^)) . (3.13) 

This serves as the definition of the pion scalar radius {r'^)s and of the coefficient Cg. 
Expanding the integral functions in Eq. (3.6) we obtain 

13 



2\7r 
S 



+xl 



-^2 (,-^ + 6^4-6L 
1 / 23 14'i 
iV V 192 ^ 2 3 ^36 



L + 



m 1 



108 



241^11 + 12{liy + Slki + nk2 - Qk4 + 6r^2 



+ 0{xl) , 



(3.14) 



and an analogous formula for 



X2 



+ 



1 19 
iVT20 
1655 1 



5 



83 



N V1152 15 



r 



15 ^ 12 ^ 



23,^ 6041 
30 ^ 2160 



1296 

3.3 Vector Form Factor 



Ki k2 

12 12 



+ fg^ 



+ 0{xl) . 



(3.15) 



As in the previous subsection we start with the general expression for the form factor. 
In the vector case it is given by 



1+X2 



-{s-4)J{s)+s{-l 



h 
6 



18N 



-X. 



-Oixl). (3.16) 
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One should notice that at zero momentum transfer gauge invariance constrains this 
form factor to be Fv{0) = 1. 

Once more we find it instructive to spht the different contributions to the form 
factor. The polynomial part of the amplitude is 



9A^ 
+s 



— (^--11 + -^fi - —L - — 
2^2^ 12 38 



47 \ 



384 192A^. 



9Nj 



23 L 



+ 



37 



36 ' 576A^ ' 86AN'^ ' ' 
and the dispersive part of the form factor is 



(3.17) 



M2) 



As) 



hli-s' + As) + - As) + hlis - 4) + h;{-s' + As) 

6 b 6 b 



7 



97 3 

s + - 

108 4 



(is' - s + 4) + ^-K,(s) (s - i) - . (3.18) 

We can now expand the form factor for s <C 1 (stt "C 4M^) and obtain the expression 



D 



where the pion charge radius is then given by 



-611 



3A;i + ^2 - \ki + 3A;6 - 12^^, 



r 

4''6 



+ 



AT 



31 



-2Z + + 



13 



181 \ 



192 48A^. 



+ 6r" 



(3.20) 



and 



2^2 . 2 

+ 2^2 



60iV 
1 



+- 



V 540 



13 ^ 1 



8429 \ 
720 ^ 25920 7 



+ r 



V2 



(3.21) 



3.4 Compcirison with earlier work 

In addition to the previously mentioned checks, we can compare to earlier partial 
results already available in the literature. 



9 



For the case of the vector form factor we have checked that the dispersive part 
agrees up to a polynomial piece with the analytical result for the dispersive part of 
Ref. [10]. In addition, we have checked that the chiral logarithms that could be 
obtained from chiral limit arguments agree with those given in Ref. [8]. 

For the scalar form factor we agree with the earlier result for the pion mass via 
the Feynman-Hellman theorem for ^^(O). We also agree with the chiral logarithms 
that were obtained in Ref. [8] and, finally we have also checked that our expression 
for Fs{s) has the correct absorptive parts as derived in [8]. 



4. Resonance and SU{3) estimates of the 0{p^) parameters 

To estimate higher order corrections due to scalar and vector resonances in both form 
factors we follow in the remainder Ref. [17], and refer to it for the notation. The 
SU(3) contributions come through kaon and eta intermediate states. These estimates 
are of course scale dependent, resulting therefore in a scale dependent final result. 
We postpone the scale dependence study to Sect. 5. 

The contribution of a given resonance state to r[ is written as with R = 

5, V, K (scalars, vectors and kaons respectively), since the effect of higher resonances 
is expected to be small [17]. 

In this section, and the for the sake of simplicity, we use a 2 by 2 matrix notation 
with 

-ir'U\ (4.1) 

the r* being the Pauli matrices and the relevant covariant derivative as defined 
in Ref. [17]. We also use 

X± = u^u^ ±ux^u, f+^u = u\v^y + a^^)u + u{v^^-a^^)v) (4.2) 

and {A) = tr(A). 

4.1 Scalar contributions 

In the scalar form factor the scalar resonance contribution is introduced via the 
lagrangian 

>C[5(0++)] = \{V^SV^S - Mis') + ca{SuX) + c^(5x+) , (4.3) 

where S contains the triplet and the singlet scalar in the leading 1/Nf. approximation. 
This corresponds to /v^ and Q = Cd/ \/2> in the notation of Ref. [17]. 

We will use the numerical values 

Ms ^ 980 MeV , c„ 42 MeV and ca ^ 32 MeV . (4.4) 
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There are more possible terms than the ones we quote in Eq.(4.3), but there is not 
enough experimental information on the scalars to determine them. Integrating out 
the scalars leads to the 0{p^) lagrangian 

^'^2lk {'''m(x+^'x+) + cl{V'iuaU^)u,u^) + 2c^c,(i.XV'x+)} ■ (4-5) 

The O(p^) lagrangian also produced in this way is included here via the contributions 
from /[ terms. There are obviously no contributions to rp, r\j, Vyi and from Eq. 
(4.5). For the rest we obtain (using F = ^ 93.2 MeV) 

'^52 ~ 'J^^™' ^^'^ ~ '^^'^^ ~ —0.3 10 , 

s 

ris ^-^CmCd^ 0.5 10-^ (4.6) 

These results should be taken as nothing more than an order of magnitude estimate 
for the size of these constants. 

4.2 Vector contributions 

Similarly as was done in [3, 4, 5, 7], we use the formalism where the vector contri- 
bution to the chiral Lagrangian starts at 0{p^). The relevant lagrangian reads 

^[^^(1^+)] = ^(VK>«1) + /x(V;K,X-]) - ^(«^) ' (4-7) 

with Vfj, describing the vector meson and V^,, — V^K — '^uV^- The parameter fy 
can be determined from p — > e'^e' [18] and gv and from p — > tttt and K* — > Ktt 
[4] , (for the latter process we must use the extension to the three-flavour case) . This 
leads to 

gv = 0.09, = -0.03 and fv = 0.20 . (4.8) 

We can now integrate out the vector meson degree of freedom obtaining only a 
nonvanishing contribution for 



_ 2V2fJ vF^ 



rY„ = ^ ^/^f « -0.25 10-3 



^v^-'-^^ 0.26 10-, (4.9) 

where we have used My — 770 MeV. 

If instead we make use of full VMD in the three-flavour case, the form factor 
would be described by 

Ml g2 4 (Mj-Ml) 
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where the last term accounts for isospin-breaking effects. We can now use as a simple 
estimate^ - = {M]^, - M^)M^/Ml, leading to very similar numerical values 



p 

for ry^ and 



"^1 ^ -0.2 10-^ and r'y^ ~ 0.21 10"^ (4.11) 



4.3 SU{3) contributions 

The SU(3) contributions to the low energy constants r[ come from kaons and etas 
intermediate states. For the vector form factor we take the expression derived within 
SU{3) CHPT [20, 21] and expand it in terms of inverse powers of M|-. This leads 
to = and 

= ^ 0.2 10~^ . (4.12) 

19207r2M| ^ ' 

The result of Eq. (4.12) is trustable but the one for Vy^ has several additional 
contributions coming from the relation between Lg and /g, the chiral logarithms and 
the relations between Fq and F (the decay constant in the chiral limit and in the 
limit with m„ = nid — rrig respectively [22]). We neglect the latter effects since 
for n !v 770 MeV the derivative of M|- log(M|-//x^) is very small. 

In the scalar form factor case we can similarly expand the expressions of SU (3) 
CHPT as given in Ref. [20] and obtain with Mk = 495 MeV 

rf, = ^ 0.3 10"^ , 

11527r2Mi 

rfa = ^ 0.9 10"^ . (4.13) 

and contributions similar to those discussed above and we also neglect 

them here, rp and only get contributions of that type, so we set both to zero. 



5. Numerical results and comparison with experiment 
5.1 General input parameters 

We use as input parameters F^ = 93.2 MeV and = = 139.57 MeV. We 

make also use of the more commonly used k quantities as defined in Eq. (3.3). We 
will mainly use the following two sets of values 

h = -1.7, I2 = 6.1, [3 = 2.9 set I , 

h = -1.5, I2 = 4.5, h = 2.9 set II . (5.1) 

more thorough treatment of quark mass corrections as done in [19] for the masses and decay 
constants would not change any conclusions within the precision needed here. 
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The value of I3 is in both cases the one derived in [2]. Set I has the other two pa- 
rameters obtained from the absolute values of the form factors using a dispersive 
improved 0{p^) CHPT calculation for three flavours [23]. Set II corresponds to us- 
ing the 0{p^) calculation of tttt scattering and fitting the D-wme scattering lenghts. 
This set also agrees, within the errors, with other determinations of h and I2 from 
dispersive analyses of tttt scattering data, see Ref. [24], and Sect. 5.2 of Ref. [4] for 
a discussion. 

In Ref. [2] a value of U ~ 4.3 was also obtained using large Nc arguments and 
the measured value of Fk/Ft^. 

5.2 The scalcir Form Factor 









0{p^) set I 


0(p^) set II 


' i 


H (GeV) 






0.5 0.77 1.0 


0.77 









0.548 


0.016 0.017 0.025 


0.054 


-0.004 


eg (GeV-^) 





5.93 


3.41 3.79 3.55 


2.52 


0.8 


Fs{0)/m 


1 


-0.0341 


0.0081 0.0086 0.0080 


0.0009 





FJF 


1 


0.0611 


0.0061 0.0061 0.0063 


0.0075 







1 


-0.0206 


0.0024 0.0026 0.0023 


-0.0003 






Table 1: Various contributions to the scalar radius (r^)^, cj, Fs{0)/{2B), F^^/F and to 
Ml/M"^. They all use h = 4.3, r[(/n) = for the columns labelled 0{p^). The quantities 
shown here do not depend on li^2 at O(p^), and therefore are not sensitive to the use of set 
I or 11 at this order. 



Although the scalar form factor cannot be accessed experimentally, indirect ex- 
perimental information can be obtained from the data on nn/KK scattering using 
dispersion relations [25] : modulo an overall normalization one can in fact derive the 
whole energy dependence of the pion scalar form factor, since it goes to zero for 
s ^ cxD as can be proven in perturbative QCD. The analysis of [25] can be used to 
determine rather precise values of (r^)^ and c^. The results for various parametriza- 
tions of the relevant phase shifts, labelled by Ai, A2 and B. give an indication of the 
uncertainty involved in such quantities. The results are - we refer to [8] for further 
explanations - 

{r^s = 0.57 W{Ai) ; 0.61 W{A2) ; 0.60 fm^B), 

cl = 10.0 GeY-^{Ai) ; 10.9 Gq\-^{A2) ; 10.6 Gq\-^{B) . ^ ' ' 
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Following [8] we will use B as central value and the range as an indication of the 
experimental error. The values labelled Ai and A2 use the CERN-MUNICH phase 
shifts [26] and those labelled B the phase shifts of Au et al. [27]. 

As was already noticed in [8], the one-loop prediction of Gasser and Leutwyler 
for the scalar radius (using I4 — 4.3 ± 0.9), 

{tYs = ^ (6^4 - y) + 0{xl) = (0.55 ± 0.15) fm^ , (5.3) 

is in very nice agreement with the dispersive evaluation of Ref . [25] . Using the same 
value for I4 and set I for the other constants, we can evaluate the corrections to the 
leading order result in Eq. (5.3) that come from the two-loop calculation (here we 
use also /i — Mp, and rg2{Mp) — 0) 

(r2)5 = 0.548 + 0.017 = 0.565 fm^ . (5.4) 

These are rather small and go in the right direction to improve the agreement with 
the experimental information Eq. (5.2). This shows that for this quantity the con- 
vergence of the chiral expansion is quite fast, and hence that it can be used reliably 
to determine rather precisely the constant I4. The central value (r^)^ = 0.60 fm^ 
given by solution B, is exactly reproduced by I4 — 4.47 (4.29), if we use set I (set 
II) together with = Mp and rg2{Mp) — 0. The influence of the latter constant on 
the value of U is tiny, and can be neglected altogether. There is some dependence 
on the choice of the scale, as it is illustrated by a variation of /x between 1 and 0.5 
GeV: I4 varies between 4.43 and 4.47 for set I, and 4.21 and 4.35 for set II. Taking 
into account the uncertainty in the determination of (r^)^, and allowing for a range 
of values between 0.57 and 0.63 fm^, we can conclude that 

[4 = 4.4 ±0.3 , (5.5) 

after averaging over what would be obtained with either set I or II (in principle, once 
the values of Zi,2 will be better determined, the error on I4 could be reduced even 
more). We stress that, given the good convergence of the chiral expansion in this 
case, the effect of yet higher orders (beyond two loops) can be safely neglected. 

The value of the coefficient at one loop does not depend on low energy con- 
stants, but it turns out to be quite far from the experimental number 

At two loops the situation improves considerably, although the exact value of this 
coefficient now depends both on 0{p'^) and 0{p^) low energy constants. Using the 
value of I4 determined with the scalar radius, we get the following numerical values 

c5(set I) = 9.85 GeV"^ + xlr's^iMp) , 

c5(set II) = 8.59 GeV"^ + xlr^ssiK) ■ (5-7) 
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Figure 1: The real part of the normahzed scalar form factor for the two sets of parameters 
that reproduce the scalar radius and, Cg of set B, for the two-loop case {0{p^)) and the 
one-loop case 



The two-loop correction to the leading order result is therefore of the right size 
needed to bring the theoretical number very close to the dispersive determination. 
To get perfect agreement we could just tunc rg^{Mp) accordingly, and get rg^{Mp) = 
5.6 X 10"^(1.5 X 10"'') for set I (set II). This fine tuning of the C(p^) constant is not 
especially interesting, were it not for the fact that the value we get is rather close to 
what we obtained with the resonance saturation hypothesis, Eq. (4.6, 4.13). Notice 
that the value of rg^{Mp) varies by a factor of three when using set I or II, and we 
should be rather satisfied with an order of magnitude and sign agreement. 

If we simply assume naive scalar dominance for the Z[ and r[ contributions to 
the scalar form factor 



Ml 



Ml 



(5.8) 



we would have obtained r^g = 8.2 10 ^ and II — 9.0 10 ^ using a value of Ms — 
980 MeV. Notice that these values are not so far from the observed ones, being of 
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Figure 2: The imaginary part of tlic normalized scalar form factor for the two sets of 
parameters that reproduce the scalar radius, and cj, of set B for the two-loop case {0{p^)) 
and the one-loop case {0{p'^)), which is identical for both. 



the right order of magnitude^. A scalar mass of 500 MeV would have increased the 
values of if and r^g by a factor of 3.8 and 15 respectively, bringing them far out of 
the region determined from experiment. 

In Fig. 1 we have shown the real part of the scalar form factor at 0{p'^) and at 
0{p^) for set I, /X = 0.77 GeV, U = 4.472, 



r 



52 



and = 4.9 10 ^, as well as 
for set II, fi = 0.77 GeV, k = 4.29, r^2 = and r^g = 16.4 10-^ In Fig. 2 we have 
shown the imaginary part for the same approximations, which is of course identical 
for the two sets at 

5.2.1 Fs(0), F^/F and M^/M^ 

In table 1 we also show the various contributions to the value of the scalar form 
factor at zero momentum transfer compared to its lowest order value. As can be 
seen, the 0{p^) correction is quite small here. For set \, ji — 0.77 GeV and l,i = 4.4 
the value is Fs{0)/{2B) — 0.974 + xlvg^. Fs{0) is of course most sensitive to the 

value of U = corresponds to II = 6.2 10"^ at = 0.77 GeV. 
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value of Is on which we have no extra information here. For completeness we have 
also included in Table 1 the corrections to the ratios Fj^/F and M^/M'^. If we use 
the value I4 = 4.4 calculated assuming the value of derived from scalar exchange 
we obtain instead 

F M'^ 

-y = 1.069 ± 0.004 and = 0.982 . (5.9) 

The error is determined by looking at the variation in table 1. We do not quote an 
error on M^/M^ since we have no improved information on I3 here. 

5.2.2 Modified Omnes representation. 

The unitarity condition which must be obeyed by the scalar form factor is satisfied 
by the following explicit representation, which is due to Omnes: 

Fs{s) = P(s)e^°(") , (5.10) 

where 

TT 74M2 S' S' — S — le 

and -P(s) is a polynomial which, in the case of the scalar form factor, can be taken 
to be a constant. In principle, if one would know the phase and inelasticity of tttt 
for 7 = 0, S'-wave, one would know also the scalar form factor up to a constant. 
Since CHPT provides a representation for the phase 5q, one could use Eq. (5.10) 
to exponentiate the result obtained at any given order of the chiral expansion. As 
discussed in [8], however, there are problems in carrying through this procedure, 
particularly because of the bad high-energy behaviour of the chiral representation 
for the phase shifts. 

The authors of [8] have instead proposed what they called the "Modified Omnes 
representation" (MOR), which is defined in the following manner. First they defined 
the reduced form factor r^(s), as 

pCHPT ^ gAA(.)pA(^^ , (5.12) 

with 

TT J AMI S' S' — S — le 

The reduced form factor has the following analytic properties: 

1. it is analytic in the complex plane, except for a cut along the positive real axis 
starting at s = 16M^; 

2. it satisfies the dispersion relation 

T\s) = l + tr dj^MlW. . (5.14) 
TT ii6M2 s' s' — s — ie 
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3. in the region 16M^ < s < Mj^ Imr^(s) only gets contributions from many 
particle intermediate states and is hence small. 

Given these properties it is easy to show that at low energy r^(s) can be well ap- 
proximated by a polynomial. Inserting a polynomial of a given order for r^(s) on 
the right-hand side of Eq. (5.12), one obtains the Modified Omnes representation. 

Gasscr and MeiBner have then compared the MOR which they obtained using a 
hnear polynomial for and the phase shift to one loop. Since the CHPT phase is 
now known to two loops, we can check here what kind of improvements this yields 
for the MOR. Besides this we also use a quadratic polynomial for F^, and fix its 
coefficients such that we reproduce the Taylor expansion at s = as given by the 
exact solution. The results are shown in Fig. 3, where we compare the MOR with 
one- or two-loop phase shifts to the exact solution. As can be seen from the figure, 

Re(F^(s)/I^i;0)) 




Figure 3: Comparison of the Modified Omnes representation (MOR) for the scalar form 
factor with the exact solution B. For the MOR we have three curves: two are made with 
set 1 for the LEC's, and one- or two-loop pion phase shift. The third one (identified by 
the word "fit") has been produced with the phase shift at two loops, and by choosing ad 
hoc values for the constants h = —1.5, h = 5.0, so as to reproduce as closely as possible 
the exact solution. 

the qualitative features of this representation do not change whether one uses the 
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one- or two-loop expression for the phase. There is of course a quantitative change 
which, however, is not very large up to 700 MeV. On this basis one could argue that 
with this representation and the two-loop phase as input, one can get a rather good 
description of the form factor up to ~ 600-700 MeV. The curves relative to set I, 
that are shown in the figure, seem to contradict this statement as their comparison 
to the exact form factor is not very good already immediately after threshold. On 
the other hand the situation improves drastically if we use set II for the LEC's and 
the two-loop phase. As can be seen from Fig. 3, an impressive agreement up to 700 
MeV can be obtained by changing I2 by half a unit with respect to that of set II, i.e. 
I2 = 5.0. We do not want to emphasize too much this agreement, since we do not 
discuss in detail the uncertainties involved. On the other hand we find it interesting 
that also the scalar form factor analyzed in this manner seems to indicate the need 
for a lower value of I2 than what was given by the K14 analysis [23]. In addition, 
this example just shows very clearly how one can use unitarity to improve the chiral 
representation, and push it somewhat beyond its typical limits of validity. 

5.3 Comparison with data for Fy 
5.3.1 The data and previous analyses 

The pion form factor has been measured both in the timelike, s > 0, and in the 
spacelike, s < 0, region. In the spacelike region there are two experiments with a 
reasonably large data set, Dally et al. [28] and NA7 at CERN [29]. The latter set 
is an accurate measurement of vre elastic scattering and dominates all fits. It agrees 
with [28] in the overlap region, where it has significantly smaller errors. 

In the timelike region there are more experiments but none of them has a large 
and accurate data set in the region relevant here. A review of the data before the 
recent inclusion of r-decays can be found in Ref. [30]. The data are obtained in 
three ways: r-decays to TnTUj- [31], e'^e~ — > tt'^tt' in electron positron colliders 
[30, 32, 33, 34] and e+e" tt+tt" measured in NA7 [35]. 

One value of the pion charge radius was obtained in [29] using a pole fit leaving 
the normalization free, leading to the result^ 

(r^)^ = 0.431 ± 0.010 fm^ n = 0.995 ± 0.002 . (5.15) 

They also used the parametrization of the vector form factor of Ref. [36] , which is a 
Fade approximation to the Omnes formula using the Si tttt phase shift, and obtained 

{r% = 0.439 ± 0.008 fm^ . (5.16) 

In [21] the same data were fitted to the one-loop CHFT formula obtained in [22] for 
the three-flavour case. There, the values (r^)y = 0.392 fm^ and (r^)y = 0.366 fm^ 

^See below for the definition of n 
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were obtained from a fit with normalization one and free respectively. The were 
not as good as for the pole fits and in the latter case the normalization ended up 
outside the error band given in [29]. The main cause was that the one-loop chiral 
formula did not satisfactorily describe the higher |s| data. 

In Table 2 we have shown the results of a fit to various sets of data (as specified 
in the table and explained below in subsection 5.3.3), using Eq. (3.19), and also a 
pole formula like 

Ms) = j^^^^ils^. (5.17) 

Notice that in this case 

Only for the data of [29] we have taken a normalization uncertainty into account. 
For none of the other experiments is the systematic error even close to the statistical 
error, with the possible exception of [35]. E.g. the CMD data of [30] have less than 
a 2% systematic uncertainty at every point. So in all fits below and in the next 
subsubsection we multiply Fy by n for fitting data from [29], and by one for all the 
other data sets. In the experiment n = 1.0000 ± 0.0045. 

5.3.2 Comparison of (r^)^ and Cy with CHPT at two loops 

The charge radius of the pion has been used by Gasser and Leutwyler to determine 
the low energy constant Iq with the result: Iq — 16.5 ± 0.9, that reproduces = 
0.439 ±0.03 fml Since we do not have other sources of information on Iq, CHPT to 
two loops can only be used here to refine the determination of this constant. It is 
instructive to rewrite the two-loop correction in the following form 




(5.19) 



where we have defined a new constant 
^6 = ^6 + 6x2 

which differs from Iq by a scale-independent quantity. From this expression it is 
clear that (besides the last piece in Eq. (5.19), which is a tiny effect) the two- 
loop correction to the charge radius consists of two main contributions. Part of the 
correction is due to the renormalization of F ^ F^^ in the leading term, and produces 
the factor 1 + 2x2/Nl4^, which numerically represents a 12% correction. The other 
part is a pure two-loop effect, and shifts Iq into Iq. Numerically, this last effect is as 
follows (using set I) 

4 - k = -0.91 + 6Nx2r'y^iMp) = -1.44 , (5.21) 
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(5.20) 



Param. 


data set 


xVdof n {r^y (fm2) (GeV-^) 


Polynom 
Eq.(3.19) 


[29] 
[29] (cut) 
Spacelike 
Timelike 

All 


42.4/42 0.995 ±0.002 0.420 ± 0.019 2.4 ± 0.5 
17.8/22 1.000 ±0.005 0.478 ±0.056 5.1 ±2.7 
50.5/55 0.996 ±0.002 0.429 ± 0.016 2.6 ± 0.4 
23.1/20 irrelevant 0.189 ±0.098 10.6 ± 2.3 
87.7/77 0.999 0.459 ± 0.009 3.5 ± 0.2 


Pole 
Eq. (5.17) 


[29] 
[29] (cut) 
Spacelike 
Timelike 

All 


41.2/42 0.997 ±0.003 0.442 ± 0.022 3.85 ± 0.68 
17.7/22 1.000 ±0.005 0.488 ±0.059 6.18 ±3.1 
48.9/55 0.997 ±0.002 0.447 ±0.018 4.0 ± 0.57 
23.2/20 irrelevant 0.191 ± 0.105 10.4 ± 2.8 
76.7/77 0.995 ±0.002 0.427 ±0.007 3.36 ±0.22 


CHPT 
Eq. (3.16) 


[29] 
[29] (cut) 
Spacelike 
Timelike 

All 


41.8/42 0.996 ±0.002 0.431 ± 0.019 3.20 ± 0.51 
17.8/22 1.000 ±0.005 0.482 ± 0.056 5.59 ± 2.8 
49.7/55 0.996 ±0.002 0.438 ± 0.016 3.35 ± 0.44 
22.9/20 irrelevant 0.134 ±0.098 11.4 ±2.3 
84.2/77 0.998 ±0.002 0.448 ±0.009 3.68 ±0.24 


Eq. (3.16) 


All 
All but [35] 


80.7/76 0.996 ±0.002 0.437 ±0.011 3.84 ± 0.25 
58.2/72 0.997 ±0.002 0.453 ±0.011 4.45 ± 0.28 



Table 2: Various fits to the pion form factor data using the simple parametrizations Eq. 
(3.19), for the first five rows, the pole formula Eq. (5.17) for the second five rows, and 
the full two-loop CHPT expression, Eq. (3.16), for the remaining seven. The last two 
fits include an extra parameter of the form d^s^ added to the two-loop CHPT expression, 
which the fit determines to be: dy = 3.0 ± l.GGeV"^ when all data are fitted, and dy = 
4.1 ± 1.6GeV-6 when all data but those of Ref. [35] are fitted. In the second of the five 
data sets used, we applied a cut and used only those data satisfying sf—s < 300 MeV. The 
errors are those that change the by one. All have a free normalization for the data of 
[29]. See text for details. 

where the numerical value after the last equal sign has been obtained inserting for 
ryi{Mp), our VMD estimate Eq. (4.9). (The use of set II shifts the above values by 
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+0.17). Modulo the uncertainty coming from the contribution of the 0{p^) LEG, 
this effect is of the order of —10%. 

The two main effects of the two-loop correction, both around ten percent level, 
contribute with opposite sign, and tend to cancel each other, resulting in a rather 
small correction to the radius 

{r%{set I) = 0.440 + 0.032 + 60.3 • rlr^{Mp) = 0.457 fm^ , 
{r%{set II) = 0.440 + 0.037 + 60.3 •r;;i(M^) = 0.462 fm^ , (5.22) 

where the final value has been obtained using our estimate for the constant ry^(Mp), 
Eq. (4.9). This result shows that also here, the chiral expansion is converging 
rapidly, and that the determination of le through the charge radius is in principle 
rather rehable. On the other hand, in this case the new LEG at 0{p^) plays a more 
important role than the corresponding one for the scalar radius case, where the size 
of the new LEG suggested by resonance saturation gave a negligible contribution. 
It is clear that lacking an independent source of information on r5)]^(Afp), wc can 
use the data on the charge radius only to determine the constant Iq. Any attempt 
to determine only Zg will depend on estimates and/or assumptions on the value of 
ry^{Mp), at least at the level of a ±10% uncertainty. 

In the case of Cy, the contribution of the O(p^) LEG rY2i-^p) ^ven more 
important 

= 0.62 + 1.96 + 1.3 X 10%2{Mp) = 5.4 GeV"^ , (5.23) 

where the first factor refers to the one-loop contribution, and the second to the two- 
loop one with ry2{^p) — 0) evaluated using the old value Iq — 16.5 and set I, and 
where the final number has been obtained using our VMD estimate Eq. (4.9). The 
coefficient Cy is therefore mainly sensitive to the value of ry2i-^p)^ ^'^'^ i'^ principle 
can be used to determine this 0{p^) LEG with rather small uncertainties (modulo 
higher-order contributions) . 

5.3.3 Fit to the data with CHPT at two loops 

Besides comparing the vector form factor GHPT formulae with the "experimental 
values" of the Taylor expansion coefficients at s = 0, we can attempt a more am- 
bitious use of our two-loop results, and directly fit the data. The reason for this is 
twofold: first, there are abundant and accurate data precisely in the region of energy 
where CHPT can be trusted more. Second, in that region, CHPT is certainly less 
model dependent than other parametrizations used to fit the data (see above), that 
make various kind of assumptions. The only assumption made in using a GHPT ex- 
pression is in the truncation of the expansion to a given order - in the present case, 
the only theoretical bias comes from neglecting contributions of three-loop order in 
the chiral expansion. As we will see, we can easily estimate such higher order ef- 
fects and therefore obtain a reliable model-independent value for the two parameters 
(r^)^ and c^. 
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Our fit has been made using tlie expression: 

Fvis)^l + Ur^)ls + cls^ + fvis) + Ois') , (5.24) 

where /v'(s) can be easily obtained from Eqs. (3.16,3.18). As free parameters in the 
fit we have used {r^)yS and Cy. It is clear that the advantage of CHPT is that it 
provides a way to calculate explicitly the function fv{s), whereas in all other cases 
one can only make guesses about its value. We remark that although the exact 
numerical value of this function depends on the LEC's, their effect is rather small, 
and the uncertainty in the knowledge of their value can be neglected altogether. 
The results of our fits are presented in Table 2, where we have fitted, as before, five 
different data sets, as specified in the Table. As can be seen there, the results of 
the CHPT fits are rather close to both those obtained with the simple polynomial 
parametrization, Eq. (3.19), and to those obtained with the pole formula, Eq. (5.17). 
This shows that the assumptions made in the two different paramctrizations, Eq. 
(3.19) and Eq. (5.17), are reasonable, and can in fact be partly justified on the basis 
of CHPT. 

To estimate yet higher orders in CHPT we have chosen to introduce an extra 
term in the polynomial part of Eq. (5.24), of the form dyS^, and to fit this new 
parameter from the data. The value we find is dy = 3.0 ± 1.6 GeV~^ when we fit all 
the data, and dy — 4.1 ±1.6 GeV~^ when we fit all the data except those of Ref. [35]. 
The changes in the values of the charge radius, and of rather small, as can 

be seen in Table 2. Moreover, the improvement in the is visible but not dramatic. 
This shows that there is no need for a large contribution from higher chiral orders, 
and that we can confidently use this extra parameter dy to estimate their effect. All 
in all, we conclude that the best values for the charge radius and Cy are given by 

{r^)l = (0.437 ±0.016) fm^ , 

4 - (3.85 ± 0.60) GeV"^ , (5.25) 

where the error also takes into account the theoretical uncertainty (i.e. we have 
added in quadrature the statistical error coming from the fit, and the difference in 
the central values between the fits with and without the cubic term in the polynomial 
- for the error on Cy, however, see below). 

In Fig. 4 we have plotted the available experimental data in the spacelike region 
together with the results of the two fits corresponding to the last two rows of Table 
2. The curves corresponding to the same two fits have been plotted in the timelike 
region, together with the experimental data, in Fig. 5. Notice that the improves 
significantly if the timehke NA7 data [35] is not considered. Our beheve is that in 
the latter experiment the systematic errors in the timelike region are underestimated. 
For this reason we have made a fit including all data but those of Ref. [35] , to show 
how much the central values would be shifted. The shift in the charge radius would 
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still be within the error bars, while that for Cy not: in view of this we have enlarged 
the error bar of this quantity accordingly (see Eq. (5.25)). 
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Figure 4: The space like data for the vector form factor and the curves corresponding 
to the fits in the last two rows of Table 2, i.e. using an expression like Eq. (5.24) with the 
addition of a cubic term in the polynomial. 

We can now use the above charge radius and Cy determinations to extract the 
values of the CHPT LEC's /g and ry2- have already stressed, one can unam- 

biguously determine from the charge radius only the constant le 

Iq = 14.6 ± 0.5 , (5.26) 

and if we want to translate this into a value for Iq, we have to make an estimate on 
the value of r^^. Using our VMD estimate, Eq. (4.9), we get 

k = 16.0 ± 0.5 ± 0.7 , (5.27) 

where the last error is purely theoretical, and takes into account the uncertainty due 
to the ryi estimate (which we assume to be ±100%), and the uncertainty in our 
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knowledge of /i 2, which also enter /g. Compared to the one-loop determination of 
the same constant, made in [2], we have found a very similar central value (due to 
the smallness of the two-loop correction we have calculated) , but we are able now to 
make a more reliable error estimate. As we have discussed, there is no way to reduce 
the error indicated here, which is mainly theoretical. Finally, from the Cy value, Eq. 
(5.25), we get 

r^2(Mp) = (1.6 ± 0.5) • 10-^ , (5.28) 
which is in reasonable agreement with our VMD estimate Eq. (4.9). 
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Figure 5: The timeUke data and two fits as indicated in the caption of Fig. 4. 



5.3.4 Modified Omnes Representation 

In principle we could have tried also for the vector form factor the use of a Modi- 
fied Omnes Representation, as we did in the scalar case: we could have even used 
that representation to fit the data. It is clear however that in the vector case the 
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polynomial part is the largest one, and that exponentiating the small dispersive con- 
tribution would not make a big effect. We have actually checked this explicitly and 
found the effect of the exponentiation to be rather small up to 700 MeV. For this 
reason we have chosen to estimate the effect of higher orders by including an extra 
term in the polynomial, as described in the previous subsection. 

5.3.5 Hadronic Contribution to the Muon [g — 2) and to a{M^) 

The process e~^e~ tt+tt" is by far the dominant part of the hadronic cross section 
at low energies. As was shown in the previous subsections, we can get a good fit 
to the data from Chiral Perturbation Theory up to about 0.5 GeV energies. We 
can therefore use our results for an improved estimate of the low energy hadronic 
vacuum-polarization contributions to the muon anomalous magnetic moment and 
q;(M|). For a recent determination and more extensive references see [38]. The 
relevant formulas are 



aMl ^ , R(s) 

AauadiMl = -— ^ Re / ds- ^ , 

6TT Jam^ s{s — M|) — le 

3 s 

R{s) = - — ^ctotl^"*^^" ~^ hadrons) . (5.29) 
The function K{s) is given in [38]. As mentioned above, at low energies we have 



(7(e+e- ^ hadrons) ^ a{e^e- tt+tt") = ^- '—^ F^{s) . (5.30) 

oS 



The total contribution to is dominated by the p region with a significant fraction 
from below 500 MeV. The contribution to the various quantities as function of the 
cutoff A^ on the integrals in Eq. (5.29) is given in Table 3 for our result for the 
CHPT form factor using the fit including all data and the d]^s^ term. In brackets we 
quote the same result but for the fit without the timehke NA7 data. The difference is 
a reasonable estimate for the error involved. This should be compared to the total 
results from Ref. [38], a^^""^ = (695±7.5) • 10"^° and Aa'^"'^(M|) = (277.8±2.6) ■ lO"! 
From the present analysis the error of coming from the region below 500 MeV is 
about 3- 10^^'^, comparable to the error on the light-by-light scattering contribution 
[39]. Once the p-mass region is better explored, more work on both the low energy 
contribution and the light-by-light scattering one will be needed. 

6. Conclusions 

In this paper we have calculated the pion scalar and vector form factor to next-to- 
next-to-leading order in Chiral Perturbation Theory and presented simple analytical 
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A (GeV) 


^QlO . ^had 


10^ • Aa(M|) 


0.32 


2.38(2.43) 


0.039(0.040) 


0.35 


7.4(7.6) 


0.13(0.14) 


0.40 


20.0(20.6) 


0.42(0.53) 


0.45 


35.7(37.2) 


0.86(0.89) 


0.50 


53.6(56.3) 


1.45(1.53) 



Table 3: Contributions of the two-pion production to a^°''^ and Aa{M2;) as a function 
of the cut-off A. 

expressions for all the relevant quantities. In addition, we have presented the known 
formulas for and using the same notation. 

We have made a careful comparison of these formulas with the data. For the 
scalar form factor this involves a comparison with the form factor derived using 
dispersion theory and chiral constraints from the tttt phase shifts as done in Ref. 
[25] . The CHPT formula fits well over the entire range of validity. Moreover, we have 
shown that by using the "modified Omncs representation" as proposed in Ref. [8] 
and which aims to resum yet higher orders by exponentiating part of the unitarity 
correction, one can improve the chiral representation, and follow quite closely the 
exact form factor up to about 700 MeV. 

For the vector form factor we have collected all available data of reasonable 
precision and performed first the standard simple fits to the data sets. Afterwards, 
we have used the CHPT formula at two loops together with a phenomenological 
higher order term to obtain a new determination of the pion charge radius and Cy-. 

{r^)l = (0.437 ±0.016) fm^ , 
4 = (3.85 ± 0.60) GeV-^ . 

The error we quote is a combination of theoretical and experimental errors, it covers 
the variation of the input parameters over the various fits and inputs done using the 
two-loop CHPT formula. 

By comparing the Taylor expansions of the measured form factors, and of their 
chiral representations, we have been able to better determine some of the LEC's that 
appear in these quantities: two of them are the O(p^) constants U and Iq, for which 
we obtained 

[4 = 4.4 ±0.3 and = 16.0 ± 0.5 ± 0.7 , (6.1) 

where the last error in Iq is purely theoretical, and where we have taken the estimated 
values of ry^ and into account in the values given. Notice that U is practically 
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free from theoretical uncertainties, as we have shown. The new value of /g together 
with the 0{p^) results quoted in [5] leads to 

[5 = 13.0 ±0.9 . (6.2) 

The other two LEC's that we have determined are 0{p^) constants, that contribute 
to the quadratic term in the polynomial of the scalar and vector form factors. We 
found 

r^3(M,) 1.5 • 10-^ , T\,,{Mp) ^ 1.6 • 10"^ , (6.3) 

with a substantial uncertainty. We find it interesting and encouraging that these 
values are rather close to the estimates we have made on the basis of the resonance 
saturation hypothesis. This result gives support to the idea that this hypothesis 
should work even at order of the chiral expansion. 
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